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ABSTRACT 



(N : a new theoretical framework for computation of the relic abundance of cold dark 

matter particles such as LSP is presented and some generic features of new results are 
discussed. The most important is a generalization of the Boltzmann equation to include 

*^ off-shell effects and its thermal average over cosmic medium. It is shown that at very low 

temperatures, much below the mass of annihilating particles, equilibrium abundance is 

Qh! 

| suppressed only by powers of temperature instead of the exponential Boltzmann factor. 

This change gives a larger relic abundance when heavy particles freeze out at these low 
temperatures. 

X 

.«?. 

1 Introduction 

The well-known problem in cosmology, the dark matter problem, requires some sort 
of new heavy stable particles which have decoupled from the rest of the universe at 
an early epoch. The conventional estimate of the relic abundance ||] of dark matter 
particles that have escaped from pair annihilation is based on the Boltzmann equation 
for the annihilating particle for which the rate is thermally averaged over the rest of 
light particles in cosmic medium. The freeze-out temperature of pair annihilation may 
roughly be estimated by equating the inverse pair annihilation rate to the Hubble rate, 
and the relic abundance is determined by the thermal number density at that freeze-out 
temperature |^], which is exponentially suppressed at low temperatures. 

t Talk given at COSM099, held at ICTP, Trieste, September 27-October 2, 1999. To appear in the 
Proceedings (World Scientific). 
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There are two possible issues as to whether this procedure gives a reliable result 
for the relic abundance. First, the Boltzmann equation, although very appealing on 
intuitive grounds, uses on-shell quantities such as the cross section and the decay rate 
(more generally S-matrix elements that relate a state at infinite past to that at infinite 
future), while quantum mechanics requires the more general transition amplitude at 
finite times, or the Green's function. The second source of concern is the use of the 
ideal gas form of the distribution function or the occupation number l/(e E ^ T — 1) at 
the freeze-out, which is very small at E r* M 3> T with M the mass of the pair 
annihilating particle. Both of these concerns are overcome and the Boltzmann approach 
is indeed justified at high enough freeze-out temperatures, as explicitly demonstrated 
in our works explained here. But at very low freeze-out temperatures, say T < Af/30, 
this conventional procedure may well be questioned. 

The key idea towards a more rigorous formulation is separation of a small system, 
in this case the system of pair annihilating particles, from the rest of a large thermal 
environment, including the annihilation product. By integrating out the environment 
part one may derive a dynamical equation for the small system and may obtain the 
quantum kinetic equation for the number density of the annihilating particle. If one 
works out the kinetic equation in a closed form and analyzes the equilibrium number 
density, this could replace the thermally averaged Boltzmann equation. The problem 
thus raised is similar to the quantum Brownian motion in thermal medium. 

The report presented here is a short summary of our recent works [[| — 1| . Before we 
go to the pair annihilation model, we explain a simpler toy model of coupled harmonic 
oscillators || , which clarifies the essential point of the temperature power law. 



2 Linear open system 

We expect that the behavior of a small system immersed in thermal environment is 
insensitive to the detailed modeling of the environment and the form of its interaction 
to the system. Since the pioneering work of Feynman- Vernon j?j and Caldeira-Leggett 
the standard model uses an infinite set of harmonic oscillators (denoted by Q(u>)) 
for the environment and a bilinear interaction with a small system (denoted by q); 

Lq = \J dio (QV)- W 2 Q 2 H) , L int = -q J diu c(cu)Q(u) . (1) 

We take for the system an unstable harmonic oscillator; the potential V(q) — \^q 2 
with loq > lu c the threshold energy. This instability condition is imposed to mimic the 
unstable particle decay and the pair annihilation process which has a threshold. 

This standard toy model is exactly solvable, both by the operator method and by 
the Feynmann- Vernon influence functional method We only quote the most essen- 
tial part of the solution for our discussion; the kinetic equation and the equilibrium 
occupation number / eq are given by 

|-rtf- W , ,„ = jyj^, (2 ) 
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H{ U ) » — < u \ , , r( W ) = ^ • (3) 

This result is obtained in the Markovian approximation which erases initial memory 
effects. The occupation number is defined here by the thermal average of the number 
operator of the system harmonic oscillator over environment variables; f(t) = (c'(i)c(i)). 
We dropped a minor (in the weak coupling limit) term of the frequency shift. The 
relaxation rate is given by T rts 7rr(o> )/ w o- 

Let us look into the equilibrium occupation / eq as a function of temperature T. 
First of all, the physical meaning of the weight function H(w) is that it is precisely the 
overlap probability of the original dynamical variable q and the eigen variable Q{w) that 
diagonalizes the total Hamiltonian. Thus, a continuous integral over u> appears for the 
equilibrium value. Indeed, this formula for / eq can also be understood by the average 
of the occupation number using the Botzmann-Gibbs ensemble of e~ Htat ' T with H to t 
the total and not the free part of the Hamiltonian. 

We consider a weak coupling region, T -c min. (uo , T) for the rest of our discussion. 
Two important cases are worth of close examination. The first is high temperature 
region in which T > lu . In this case the Bose- Einstein distributiion function l/(e"/ T — 1) 
hardly varies where the Breit-Wigner function H(u>) changes over T, and one may replace 
H{lo) by the delta function S(uj — u>o)- This gives / eq « l/(e"°/ T — 1) , which is nothing 
but the on-shell form of the occupation number. On the other hand, if the temperature 
is very low like T -C wq, the region around to ~ loq gives a minor contribution in 
the continuous integral, and the threshold region u ~ uj c becomes more important. 
Parametrizing this region by r(u) ~ c(u—u c ) a , one gets / cq ~ 0[1] x cT Q+1 /(w 2 — w 2 ) 2 . 
The salient feature of this formula is appearance of the temperature power law oc T a+1 ; 
it arises from the continuous u integral outside of the on-shell region u> ~ wq. Even at 
intermediate temperatures the sum of two limiting formulas is a good approximation; 

/c ^^^ + ° [1]x R^- (4) 

Tha basic features of this simple model can be taken over to the more complicated 
cases such as the unstable particle decay and the pair annihilation to which we now 
turn to. 



3 Pair annihilation model 

We consider a simple boson pair annihilation model. The heavy particle is denoted 
by the field <p and the light environment particle by X- The interaction Lagrangian 
density that gives rise to the annihilation process ipip — ► \X lii 

C int = -^ 2 X 2 ~^V 4 -^X 4 + SC (5) 

The last two terms (oc X v or A x ) arc introduced both for consistency of renormalization 
and for thermalization of x particles, hence the parameters are taken to satisfy \\ v \ <C 
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A 2 C 1, |A X | < 1 , but A x | 3> |A|. The renormalization counter term 5C is determined 
in the usual way. 

It is convenient for discussion of our approximation scheme to employ the influence 
functional {jj. This quantity denoted by J-[q(r) , q' (t)] results after the environment 
integration for the density matrix of the entire system. Since the density matrix involves 
both the transition amplitude and its conjugate, the path integral formula resulting from 
the environment integration has a functional dependence both on the path q(r) and its 
conjugate path q'(r). For the environment field X we may take the Gaussian part of 
the Lagrangian density since the A x self-interaction only serves for thermalization of the 
X system. With the thermal ensemble given by a Gaussian x density matrix, one can 
explicitly perform the x path integration. The result is a nonlinear quartic action; 

exp[-— / dxdy (b(x)a R (x- + i&(x)ai(x- y)X 2 (y))} , 

X 2 (x) = p 2 (x)+ip' 2 (x), £ 2 (x)=<p 2 (x)-p' 2 (x), (6) 
a(x) = a R (x)+ia I (x) = \ 2 (tr (t[ X 2 (x) X 2 (0) p<*>]) - (tr X 2 p^ ) 2 ) , (7) 

(y) 

with pp the thermal density matrix for the environment x- The last equation for the 
kernel a means that it is the real-time thermal Green's function. 

One point is important before we further go on. Since we deal with a non-equilibrium 
state for the ip field, time translational invariance is in general violated; for the correlator 

(tp(x,t 1 )tp(y,t 2 )) ^ (p(x,ti -t 2 )(p(y,0)) , (8) 

unlike the correlator in complete thermal equilibrium. Thus, the correlator has depen- 
dence both on the relative and on the central time. 

It is crucial for further development to realize that dependence on these two times 
is greatly different. Variation in relative times is governed by the inverse of the system 
frequency 1/loq and is fast, while the central time dependence is given by a much larger 
time scale, l/(A 2 Wo) since without interaction the system does not change. It is thus 
reasonable to treat the central time dependence adiabatically; one first regards the 
central time as a constant and solves the dynamics with respect to the relative time, 
and finally recovers the central time dependence at the end. 

Another technical point for simplification is the Hartree approximation. The influ- 
ence functional given above has four fields ip at two different times in the exponent. 
In a dilute system such as the pair annihilation in the expanding universe it is a good 
approximation to ignore higher orders of correlation. This is in accord with the spirit 
of the Hartree approximation in which one replaces a pair of p fields in the influence 
functional by the correlator yet to be determined. If one formally solves the model in the 
Hartree approximation, the resulting correlator becomes a functional of the correlator it- 
self, thus one arrives at a self-consistency equation for correlators. This self-consistency 
equation is not very convenient for further analysis and it is far better to derive from the 
self-consistency relation time evolution equation, namely the quantum kinetic equation. 
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A convenient quantity for the kinetic equation is a combination of Fourier transform 
of the correlator, 

a(k, t)= ( d\x - y) (<p(x) <p(y)) e^-a) , (9) 



with t = (xq + yo)/2. Since spatial homogeneity holds for the problem of our interest, 
this quantity a does not depend on the central spatial coordinate (x + y)/2. We define 
the key quantity r(k , t) = o{—k , t)/<7_ (k , t) , where <j± is even and odd parts. With 
the slow variation as to the central time one has the kinetic equation in the form S , 
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with E p = \/p 2 + M 2 , = vfef + ti 2 • In the right hand side r = r(p , t) , r 1 = r(p' , i). 
The distribution function in the right hand side fi is given by that of light \ particles 
in thermal equilibrium. The weight function H here obeys a self-consistency equation; 

H(k s r_(fc,oo) , . 

1 (fc 2 -M 2 (T)-n(fc,oo)) 2 + (7rr_(fc,oo)) 2 ' 1 ' 

r_(fc,oo) = 8 / — -H(k' ,00) r x (k + k')- 



( e /3(fco+fc'o) _ _ e -Pk> \ ' 



r x (uj,k) = -£—[Jl p^—e{io - \/k 2 + 4m 2 ) + In ■ 1 



xv ' ; 256tt 2 I V t^ 2 -fc 2 v ; /3fe i- e -/S|w-l J ' 

with /3 = 1/T the inverse temperature and u>± = ^ ± | \/l — 4m 2 /(aj 2 — fc 2 ) . We 
included the temperature dependent mass shift M 2 {T) of order A and further 0[A 2 ] 
terms of the proper self-energy II. 

The structural resemblance of this kinetic equation to the Boltzmann equation is 
obvious and in this sense this equation is its generalization. There are a few important 
differences, however. First, the basic quantity r(p ,t) is a function of 4-momentum p 
unlike the distribution function f(p,t) in the Boltzmann equation, which only depends 
on the 3-momentum. Related to this is that even the processes not allowed by the 
energy-momentum conservation on the mass shell contribute to the collision term in the 
right hand side, for example 1 «-> 3 process such as ip <-> <pxX- 

Obviously, the equilibrium solution for which the right hand side vanishes is r(p , 00) — 
\j{eP Vo — 1) . We assume that this solution is unique and no other solution exists. One 
further determines the equilibrium correlator by 

(<p(xMy)U = J ^-J- w -H(k, (X )e-^~yK (12) 
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It can be shown that to 0[X 2 ) this correlator coincides with that given by the thermal 
held theory {§. 

The weight function H or r_ obeys the self-consistency equation. We analyze this 
and the kinetic equation perturbatively. In the zero coupling limit H(p, oo) « 6(po — 
Ep)/2E p for po > 0, and our kinetic equation reduces to the usual Boltzmann equation. 
To include the next 0[X 2 ] order the weight function H(k , oo) is calculated by using the 
integral form of r_ (k , oo) in which H is replaced by the delta function. The resulting 
kinetic equation with H{k , oo) to this order and its associated equilibrium formula gives 
our basic results. 

Once the kinetic equation is solved, one can compute physical quantities using the 
function r. The ip energy density is thus calculated taking the coincident time limit of 
two-point correlators, with due consideration of renormalization of composite operators. 
Thus, 

1 1 M 2 
(H<p) = + 2^ 7 ^ 2 + ~ ( counter terms )) ( 13 ) 

- J J d Po (p 2 + E 2 p )H(p , oo)r(p, t) 

A 2 f d 3 kd 3 k' n2 t ;/2w , p p 2 

/d 3 v T 6 
j^E p T(E p ,p,t)+c\ 2 W , (14) 

with c = ggi^ ~ 1.4 x 10" 5 . (The X mass m was taken to vanish for this calcula- 
tion.) In the last step we replaced the spectral weight H by the delta function. The 
temperature power term (oc T 6 ) thus arises, as expected. The first term for the energy 
density involves the r function on the mass shell, hence is Boltzmann suppressed at low 
temperatures. 

The temperature power T e is understood as follows. Since the operators for TL V have 
mass dimension 2 or 4, terms of order T 2 and T are divergent and they are cancelled 
by counter terms. The remaining finite term is of order T 6 , as given above. Precisely 
the same temperature power term was also derived by using the technique of thermal 
field theory Q, which is hardly surprising since our equilibrium result is given by the 
Boltzmann-Gibbs ensemble. 



4 Cosmological evolution 

We find it reasonable to define the ip number density by n v = p v /M at low temper- 
atures of T <C M. The time evolution equation for the number density in the expanding 
universe is given by 

^+3Hn ip = -(av)(nl - n 2 q ) , (15) 
with H the Hubble rate, {av) w A 2 /(32vrAf 2 ) , and 

-e q «(^) 3/2 e-^ + cA 2 ^. (16) 
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The usual time-temperature relation t oc T~ 2 may be used. Time evolution calculated 
from eq. (|l5|) supports the picture of the sudden freeze-out as in the Lee- Weinberg anal- 
ysis Q. The major difference here is the new equilibrium abundance n cq . When the 
freeze-out temperature is high enough, our result justifies the usual thermally aver- 
aged Boltzmann approach. On the other hand, when the freeze-out occurs at which 
the temperature power term dominates, the final relic abundance is enhanced over the 
Boltzmann abundance. 

For numerical results I refer to our paper Q. Our new effect tends to show up 
for a larger coupling and a smaller mass. Since the new effect gives an additional 
positive contribution to the energy integral, the relic density is always enhanced over 
the conventional result. Thus the allowed parameter region in the model parameter 
space gets always shrunk by our result. 
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